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Abstract 



We generalize a recently proposed mechanism for the origin of primor- 
dial metric perturbations in inflationary models. Quantum fluctuations of 
light scalar fields during inflation give rise to super-horizon fluctuations of 
masses and reaction rates of various particles. Reheating, freeze-out, and 
matter-domination processes become inhomogeneous and generate super- 
horizon metric perturbations. We also calculate the degree of non-Gaussianity 
f n i for this new model of cosmological perturbations. The precise value of f n \ 
depends on the specific models, but \f n l\ ~few is a natural lower bound for 
our mechanisms. This is much larger than the currently assumed theoretical 
value f n \ ~ tilt < 0.05, and is thought to be observable. In a particularly 
attractive model of inhomogeneous mass-domination, the non-Gaussianity of 
perturbations generated by our mechanism is simply f n \ = 5, irrespective of 
the detailed structure of the underlying field theory. 
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I. INTRODUCTION 



During inflation the energy density is dominated by the potential energy of a slowly 
rolling scalar field, the inflaton. At the end of inflation this energy density has to be converted 
into normal particles, reheating the universe and starting the standard phase of the hot 
big bang. In [1] we suggested that if the inflaton decay rate T varied in space, density 
perturbations would be generated during reheating independently of those generated by the 
standard inflationary mechanism. If two different regions of the universe had different Ts 
then effectively the inflaton would decay into radiation first in one region and then in the 
other. During the time one region is filled with radiation while the other one is not, the 
universe expands at a different rate in each region, resulting in density perturbations when 
reheating is finished. The decay rate of the inflaton is determined by the expectation values 
of some scalar fields. If those scalar fields were light during inflation they fluctuated, leading 
to density perturbations through the proposed mechanism. Here we extend and generalize 
this model (§11). 

We also calculate the non-Gaussianity of the metric perturbations generated by our 
mechanism. The standard inflationary model predicts unobservably small non-Gaussianities, 
while our models generically predict non-Gaussianities of potentially observable magnitude 
(PI). 

II. MASS-DOMINATION MECHANISM 
A. The Mechanism 

We first discuss a generalization of our original mechanism [1]. Assume that density per- 
turbations created during inflation, as well as during reheating, are negligible. So right after 
reheating the Universe is filled with radiation of uniform energy density and temperature 
T R . Assume that the mass and the decay rate (T) of some of the created particles (call them 
ip) is set by the VEV of a scalar field 0, 
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M = A0 (1) 

where A is a coupling constant. We shall assume that the mass of during inflation is 
smaller than the Hubble parameter. Thus, inflationary fluctuations of <fi on superhorizon 
scales get imprinted into the mass of ?/>-quanta. As long as the temperature is much larger 
than M, the perturbations in mass do not contribute into the energy density. However, once 
the temperature drops below M, the mass-fluctuations become important. For the mass 
fluctuations to be imprinted into the density perturbations, it is essential that ip dominates 
(or at least becomes a significant component of) the energy density of the Universe during 
some period. That is, the decay rate of ip must satisfy 

r < h d (2) 

where Hp the Hubble parameter at the moment when tp starts dominating the energy 
density. For simplicity, we shall assume that the annihilation rate of ip is much smaller 
than M 2 /Mpi. Then ip start dominating as soon as they become non-relativistic (T ~ M, 
assuming that ip were thermalized at some early time). The corresponding Hubble parameter 
is H D ~ M 2 /Mp t . 

If (2) is satisfied, the perturbations in M and T get translated into the density pertur- 
bations. Due to the variation of M and T, the interval of -^-domination in different regions 
of the Universe will be different leading to density perturbations at the end of the process. 

We shall now derive the magnitude of the resulting density perturbations. We compare 
final radiation energies (after ^-decay) in different superhorizon regions, for the same value 
of the scale factor a. It is simplest to compare energy density in any given region to the one 
in radiation, which scales as 

Prad oc a" 4 , (3) 

and is the same in all the regions of interest. The -?/>-energy density in a region with mass 
M and decay rate T scales as radiation for all the values of a, except the interval of its 
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domination, that is, between the domination and the decay. Domination starts at a = 

(^domination when 

P Pdomination M (4) 

and ends at a = a decay , when 

P = Pdecay = T 2 M Pl (5) 

Outside this interval a decay > a > a domination , the energy density in all the domains scales 
as radiation and is independent of either M or T. However, during the domination interval 
energy of ip scales as matter 

p oc a" 3 (6) 
and becomes different in different domains. Thus, we have 



ddecay \ Pdomination 



(^domination ' Pdecay 

The final energy density stored in ip, right before the decay is 



M A T- 2 Mpf (7) 



P oc Qdecay Prad = M^Y~^M p 2/ "p rad (8) 

(^domination 



The resulting density perturbations are given by 

5p 4 SM 2 5T 



(9) 



p 3 M 3 r 

Depending on the underlying model, one can consider several possibilities. For instance, 
for M = constant, this recovers our previous result [1], where now we are interpreting T as 
the decay rate of the inflaton rather than tp. The results carry over from one scenario to the 
other because during reheating, when the oscillating inflaton dominates the energy density 
we also have p oc a~ 3 . Another interesting case for the realistic model building is T oc M, 
for which we get 

*e = -— do) 

p 3 M y ! 
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B. Constraints 



The necessary requirement for our mechanism is that an effective mass of must remain 
smaller than the Hubble parameter during inflation. However, this relation may or may 
not be violated during -^-domination period. If it is violated, then will start oscillations 
during the ip domination epoch. The calculation of the resulting density perturbations for 
oscillating is more involved and will not be discussed here. 

The simplest situation is when the effective mass of remains smaller than the Hubble 
parameter all the way until the ip particles decay. This requirement strongly constraints 
the initial value of 0, provided the 0-dependence of -0-mass is significant. Below we shall 
estimate the constraint for the case of maximal dependence, that is when the entire mass 
of -0 comes from according to (1). Generalization for the case of milder dependence is 
obvious. 

In many cases this puts a non-trivial constraint on the initial value of 0, as we shall 
now demonstrate. The crucial point is that at high temperature the non-zero density of ip- 
particles generates a large thermal mass for 0. For instance, every scalar degree of freedom 
that is in thermal equilibrium at temperature T, and is coupled to 0, generates the following 
contribution to the 0-mass 

ml ~ A 2 T 2 (11) 

This contribution is there irrespective of itself being in thermal equilibrium. This mass 
may exceed the value of the Hubble parameter during the ^-domination, unless the following 
condition is met 

AT D < H D (12) 

where To = A0 is the domination temperature. This implies that must satisfy the 
following condition. 



Notice that (13) guaranties that the condition < H will be satisfied both before and 
throughout the domination. This is easy to understand. Above the domination temperatures 
the thermal mass scales as temperature, and quickly becomes negligible relative to Hubble 
parameter, which scales as T 2 . During the domination period ^-particles are out of thermal 
equilibrium and their energy scales as non-relativistic matter, according to (6). So does the 
mass of 0, m| oc p^. Thus, and H scale in the same way, and relation < H is 
maintained throughout the domination. 

For typical 0-dependence of M and r, ~ M Pi implies that the amount of density 
perturbations created by standard inflationary scenario will not be sub-dominant. So, our 
mechanism will be the primary source of perturbations only in inflationary scenarios that 
have no fluctuating inflaton field, e.g., such as recently proposed "self-terminated inflation" 
[2]. Even in this case the level of the gravitational wave background should be significant. 

We should again stress that the above constraint in absent in the models in which M is 
independent of 0. In these class of models our source of density perturbations can dominate 
over the standard inflationary mechanism by many orders of magnitude. 

C. Models 

We shall now consider some practical implementations of our mechanism in realistic 
models. As is obvious from the above, we are interested in theories in which masses of 
dominating particles M(0) and/or their decay rate(s) r(0) are functions of the fluctuating 
flat-direction field 0. Since, in general 0-dependence of functions M(0),T(0) can be very 
different, it is hard to talk about more than an order of magnitude predictions. However, 
there is a sub-class of theories in which all the mass-scales (at least during the epoch of 
interest) are set by alone. In such a case M(0) oc 0, T(0) oc 0, and density perturbations 
via our mechanism acquire an especially simple form, practically independent of either the 
coupling constants or the field content. Interestingly, the Standard Model, as well as its 
minimal supersymmetric extension, in which is identified with a flat direction field, fall 
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within this category. 

Before proceeding we have to stress that there can be corrections to the masses and decay 
rates from other existing scales in the theory. For instance, non-perturbative gravity effects 
can contribute into the decay rates of the Standard Model particles. If such corrections are 
there in the first place, their relative value will depend on many factors, such as existence of 
discrete gauge symmetries, masses of decaying particles etc. It is not hard within a concrete 
model to make this corrections sub-dominant. This is beyond the scope of the present work, 
in which we shall deal only with the simplest minimal case. 

We shall now apply our mechanism to Standard Model, and to its minimal supersymmet- 
ric extension. We assume that (1) during inflation at least some of the MSSM flat direction 
fields (call them 0) have masses <C H; and (2) after the inflaton decay, the Universe is 
left with a thermal gas of SM particles (and their superpartners) of uniform energy density. 
Under these assumptions we will demonstrate that the amount of the density perturbations 
generated through the cooling process, irrespective of the coupling constants or the nature 
of flat direction, is given by 

t = It- (14) 

P 3 

provided some of the particles dominate the Universe for a short period. 

To show this, let us first ignore superpartners and only consider Standard Model particles. 
Let be a flat direction corresponding to the Standard Model Higgs field. That is, we shall 
assume that the value of the Higgs mass during inflation satisfies our requirements, and 
that the Universe is reheated by producing some of the standard model particles. When 
the Universe cools down, the heavy particles decouple and decay into the lighter ones. The 
heavy states that have small annihilation and decay rates can dominate the Universe. Such 
can only be fermions (quarks and leptons), since the gauge bosons decay through the order 
one gauge couplings and decay before domination. 

For a given species to dominate it must have frozen out but still have not had enough 
time to decay. Thus both its annihilation and decay rates must be less than H. This 



situation can always be achieved for some species for the appropriate initial value of 0. 
Assume for defmiteness that ~ M P , and, of course, we shall assume that T R <C M P) but 
big enough that some of the unstable fermions are produced. We shall now discuss under 
what conditions some of the species can undergo a brief interval of domination. It is useful 
to consider electrically charged fermions and the neutral ones separately. 

Charged fermions can annihilate into photons so their freeze-out abundance will be much 
smaller than M 3 unless their mass is of order M ~ a 2 EM M P i. Their mass however has 
to be smaller than the reheating temperature at the end of inflation for them to be in 
thermal equilibrium to start with, so M < T r h ^5 \J H in fMpL ^5 10~ 2 Mpl- This case is only 
marginally possible so we will concentrate on uncharged fermions, ie. neutrinos. 

The story with neutrinos is very different, as we shall now discuss. In order to understand 
how neutrinos can dominate in very early Universe, we shall specify the origin of their 
mass. For defmiteness, we shall assume the standard "See-Saw" mechanism [3], in which 
the neutrino masses are generated by mixing to the heavy gauge-singlet fermions (right- 
handed neutrinos). The relevant terms in the Lagrangian are 

K<t>viMR + M R v R v R + ... (15) 

where M R is the Majorana mass of the gauge-singlet fermion A is an Yukawa coupling 
constant, and generation indicies are suppressed. 
In today's Universe 

\ v (f> < Mr, (16) 

and heavy neutrinos can be integrated out. As a result of this integration the light neutrinos 
acquire small Majorana masses given by 

However, in the epoch of our interest ~ M P i, and the condition (16) can be violated. 
Depending whether this is the case, the light neutrino will either continue to be a Majorana 
particle with mass (17), or will effectively become a Dirac particle of mass 



My ~ \ v (f> (18) 

Neutrino annihilation rate into the fermions happens through the exchange of heavy Z, W 
bosons, with masses ~ Mpi, and the annihilation into the Higgs particles is suppressed by 
\i/M%. So it is safe to assume that their abundance froze out early on. 

The decay rate of neutrinos requires some attention. In the early Universe some of 
the neutrinos that are stable today, could have been unstable. For instance, since the 
electron and quark masses scale linearly with as opposed to neutrino masses that scale 
quadratically, the electron neutrino u e can become heavier and decay into the electron and 
the quark- ant i- quark pair 

v e -> e~ + u + d (19) 

In general, the Standard Model fermions decay through the exchange of VU-bosons and 
the decay rate is: 

T = /M 5 /0 4 (20) 

Where / is a small constant that depends on mixing angles, and is different for different 
fermions. The crucial point is that because M is set by 0, T is linear in (for neutrino this 
requires \ u (j) > M R ) 

r = /a 5 (21) 

As a result of this, irrespective which fermion happens to dominate, the imprint of density 
perturbation is universally given by (14). It is easy to check that for ~ Mpi ~ 10 4 H, 
the domination condition is always satisfied by at least some SM fermions. For large values 
of and H, our mechanism can operate several times and each time imprint the density 
perturbations. 

Interestingly, going to MSSM does not modify this general result, irrespective which flat 
direction develops large VEV during inflation. This can be understood in the following way. 
Let be some MSSM flat direction. breaks gauge symmetry and gives masses oc to 
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some gauge fields and fermions (and their superpartners). Since, for the needed magnitude 
of perturbations the flat direction VEV must be >> H, the supersymmetry breaking effects 
can be ignored. So we can neglect the mass splittings between the superpartners. Thus, 
irrespective of the particular flat direction, the ^-dependence of T will be the same. This is 
obvious since the masses of particles are set by 0, and heavy particles decay into the light 
ones through the exchange of the massive gauge bosons. 1 This leads us to a conclusion 
that the density perturbations generated via our mechanism is independent of the detailed 
structure of the couplings as well as the nature of dominating flat directions, and is always 
given by equation (14). 

III. NON-GAUSSIANITIES 

The WMAP satellite has provided bounds for the degree of non-Gaussianity of the pri- 
mordial cosmological perturbations. Assuming that the super-horizon gravitational potential 
fluctuations (during the matter-dominated era) are of the form 

® = g + fnig 2 , (22) 

where g is Gaussian, the WMAP results are —58 < f nl < 134 at 95% confidence [4]. The 
Sloan Digital Sky Survey should provide similar accuracy [5]. 

The standard one-field slow-roll inflation predicts the degree of non-Gaussianity which 
corresponds to ~ tilt of the perturbation spectrum [6]. The power spectrum measured 
by WMAP is consistent with scale invariant, however when WMAP data is combined with 
other probes of large scale structure small tilts maybe preferred [7]. It is fair to say that the 
largest tilts still allowed by the data are of order \n — 1| < 0.05. Thus the current theoretical 
expectation is that primordial non-Gaussianity is unobservably small. We will show that in 
our model non-Gaussianities are much larger. 

1 The massless exchange cannot lead to a particle decay, due to the fact that the massless gauge 
bosons couple through unbroken generators, which commute with the Hamiltonian. 
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We will consider non-Gaussianities assuming that the fluctuations of the light field 
which is responsible for fluctuations of the coupling constants and/or masses are purely 
Gaussian. The case of non-Gaussianities of will be discussed elsewhere. 

To talk about non-Gaussianities, one needs to define a gauge-invariant quantity to 
quadratic order. We will follow [6,9,10], and use a gauge invariant variable ( that remains 
constant outside the horizon. It is defined as follows: is proportional to the local scale 
factor measured on uniform local Hubble parameter hyper surf aces. To linear order, ( is pro- 
portional to the gravitational potential <3>. During matter domination ( = —(5/3)$, during 
radiation domination ( = — (3/2)$ [11]. 

In our original scenario [1], the energy density fluctuations on hypersurfaces of constant 
scale factor a are 

p oc r~ 2/3 . (23) 

Since p oc a~ 4 , this gives 

C = -(l/6)lo g r. (24) 

The inflaton reheating rate T = A 2 m^, where is the inflaton mass at its minimum, 
and A is the coupling constant of the inflaton decay. We have to consider two case. First 
assume that fluctuations of A are negligible, and oc dominates the fluctuations of V. 
Then, up to second order, 

< = 4( f *-5*j)- < 25 > 

where 6$ = ((f)/ < > — 1) is assumed to be Gaussian. The definition of /„; used by [8] 
corresponds to [6] 

C = g- (3/5)f nl g 2 : (26) 

And gives 

fm = -5 (27) 
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for this scenario. If, on the other hand, A oc dominates the fluctuations of the reheating 
rate, one obtains a smaller non-Gaussianity, 

fni = -5/2. (28) 

This shows the general trend: less efficient mechanisms of translation of the 0-fluctuation into 
the metric fluctuation give greater non-Gaussianities (assuming the inefficient mechanism is 
still the dominant one). Our assumptions that m oc <fi and A oc <fi can also be generalized to 
include non-linear terms, leading different values of 

The new mass-domination scenario described in §11 is of particular interest because of 
its high universality, based on being the only mass scale of the theory. From p oc 2 ^ 3 , we 
get a unique answer independent of the Yukawa-couplings: 

fni = 5. (29) 



IV. CONCLUSIONS 

We proposed a new scenario for generation of metric perturbations after inflation, 
whereby superhorizon fluctuations of the light fields generated during inflation are translated 
into metric fluctuations during reheating, mass domination, or freeze-out. 

We calculated non-Gaussianities for this model. Our result is f n i = 5 for the mass 
domination mechanism. Other scenarios give different values, and non-Gaussianities of 
the light fields <fi also lead to partial compensation or enhancement of the intrinsic non- 
Gaussianities of our mechanism, but an observable value \f n {\ ~few is a natural lower bound. 
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